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Summary.  The  problem  of  reconstructing  the  density  o*  a vibrating  string  given 
the  first  N eigenfrequencies  for  two  vibrating  configurations  admits  an  infinite 
number  of  solutions.  Among  aii  such  strings  compatible  with  the  truncated  data 
set,  we  define  the  ideal  string  to  be  that  string  for  which  a weighted  average 
of  the  density  is  minimum.  We  prove  that  this  ideal  string  must  have  a finite 
number  of  degrees  of  freedom  and  hence,  that  it  is  made  up  by  a finite  number 
of  concentrated  point  masses.  By  specializing  the  optimality  criterion,  we  can 
also  show  that  the  Krein  string  is  an  ideal  string. 
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I'  i 

!'  : 

ii  1.  Introduct  ion.  " 

' j 

This  paper  is  concerned  with  the  solution  of  inverse  eigenvaiue  problems  j 

for  which  the  data  sets  are  insufficient  for  insuring  unique  solutions.  This  | 

I 

situation  is  typical  of  most  inverse  eigenvalue  problems  arising  in  Geophysics.  | 

A very  appealing  approach  to  certain  inverse  problems  with  partial  data  j 

sets  was  pioneered  by  Parker  (197^»  1975).  His  original  work  dealt  with  inverse 
problems  associated  with  the  reconstruction  of  a buried  body  embedded  in  a known  I 

matrix  from  surface  gravity  anomaly  data.  Confronted  with  the  inherent  nonunique-  I 

ness  of  the  problem,  he  decided  to  construct  that  particular  body  which  was  (i)  | 

I 

capable  of  explaining  the  data  and  (ii)  had  the  least  possible  maximum  density.  ' 

I 

Parker  called  this  unique  body  the  ideal  body.  Thus,  out  of  a set  of  many  possible  ’ 

solutions,  he  singled  out  a particular  one  by  means  of  an  extremum  criterion.  This  < 

procedure  yielded  rigorous  bounds  for  the  entire  class  of  solutions.  « 

A similar  approach  has  also  proved  very  successful  is  studying  certain 
turbulent  flows  (Howard  1963,  1972,  Busse  1969).  There  the  problem  stems  from  i 

the  fact  that  the  Navier-Stokes  equations  admit  a multiplicity  of  solutions.  ' 

Out  of  this  class  of  solutions,  a distinguished  one  is  selected  by  means  of  an 
optimality  criterion.  This  optimum  solution  can  then  yield  general  results,  e.g. 
bounds  on  all  solutions. 

The  present  paper  attempts  to  approach  the  simplest  inverse  normal  modes 
problem,  namely  that  for  a vibrating  string,  from  the  same  point  of  view.  The 
data,  which  consist  of  truncated  frequency  spectra,  cannot  guarantee  a unique 
solution  to  the  inverse  problem.  Among  the  set  of  strings  which  have  the  same 
given  eigenf requencies,  we  shall  select  an  ideal  string.  As  a selection  criterion, 
we  shall  minimize  a weighted  average  of  the  density.  We  shall  prove  that  this 
Ideal  string  has  a finite  number  of  degrees  of  freedom  and  hence,  is  made  up  of 
concentrated  point  masses  joined  by  v/eightless  threads. 
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The  particular  features  of  the  ideal  string  can  be  found  by  solving  a 
convex  programming  problem.  Sabatier  (1977a, b)  was  the  first  to  point  out  the 
relationship  between  Parker's  ideal  body  theory  and  the  classical  theory  of 
linear  programming.  A similar  relationship  exists  here.  For  a special  ideality 
criterion,  the  solution  to  this  programming  problem  can  be  trivially  obtained. 
This  special  case  is  discussed  at  the  end  of  the  paper. 
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2.  The  extremum  criterion  for  the  ideal  string. 


The  vibrating  system  of  interest  consists  of  a string  of  (dimensionless) 
density  p{x),  of  unit  length  and  taut  by  a unit  tension.  Two  vibrating  configu- 
rations will  be  used,  differing  from  each  other  by  the  fastening  at  the  left  end. 
The  corresponding  eigenvalue  problems  will  be  written  thus: 


u''  + X p u ■ 0 , I 

n n n 


u (0)  cosot  - u*(0)  sina  ■ 0 
n n 

u (1)  COSY  + u’(l)  siny  ■ 0 
n n 


V (2.1) 


and 

I I 

V + y p V - 0 , 

n ^n  n ’ 

V (0)  cos6  - v’ (0)  sin6  * 0 

n n 

V ( 1 ) cosy  + u * ( 1 ) s i ny  ■ 0 

n n 


(2.2) 


ot,6,y  are  parameters  such  that 
0 ^ a < e < tt/2. 

As  a result  the  system  represented  by  (2.1)  is  stiffen  than  that  associated  with 
(2.2)  and  the  eigenfrequencies  interlace  as  follows: 


U|  < Xj  < • • 

• y < X < • • • 

n n 

(2.3) 

preferable  to  write 

(2.1)  and  (2.2)  as  integral  equations,  viz. 

f p(C)  G(x,a  u^(C)dC 

u (x)  - X 

(2.4) 

n n 

JQ 

n n 

f’  p(c)  r(x,c)  v^(c)dc 

jf) 

(2.5) 

where 


- It  - 


cosacosY  + s i n (a+^ 


(2.6) 


{(l-C)cosY  + sinY}{xcosa  + sina} 

■ COSCCOSY  ! sTn-Sn)  ’'■'  <2.6) 

{(l-x)cosY  + sinY}{Ccosa  + sina} 

X > C 

and  r(x,C)  is  obtained  by  replacing  a by  6 in  the  above  formula. 

Given  p(x),  (2.U)  and  (2.5)  can  be  looked  upon  as  recipies  for  calculating 

the  corresponding  eigenvalues  and  y^.  For  brevity,  we  denote  the  result  of 

these  calculations  by  the  notations  X [p]  and  y„[p]  which  emphasize  that  X and 

n n n 

are  functionals  of  p(x). 

It  is  well  known  (Borg  19^6),  that  the  complete  X and  y spectra  are  necessary 

(and  sufficient)  for  the  unique  determination  of  p(x).  Therefore,  if  we  are  only 

given  the  truncated  spectra  {X  },  and  {y  },  we  cannot  infer  p(x)  uniquely:  there 

n I n I 

are  infinitely  many  strings  with  the  first  N natural  frequencies.  Let  us  denote 
the  set  of  these  strings  by  i.e. 

Rjj  - {p(x);  p(x)>0,  X^[p]  - X^,  y^[p]  » y^,  n-1 , 2,  . . . N}  (2.7) 


Next,  we  define  a weighted  average  of  p(x),  viz. 

M[p]  - 1 f(x)  p(x)  dx  , 

^0 


(2.8) 


where  f(x)  Is  a positive,  continuously  differentiable  function  of  our  own  choosing. 
The  only  other  condition  which  we  shall  place  on  f(x)  is  that: 

f(l)  COSY  + f’(l)  slnY  - 1 . (2.9) 

For  instance,  if  y * 0,  we  can  take  f(x)  to  be  equal  to  x , in  which  case 
H[p]  would  correspond  to  the  k-th  moment  of  the  density  distribution. 

We  can  now  state  the  extremum  criterion  which  will  define  the  ideal 
string.  From  among  all  the  strings  in  Ry,  the  Ideal  string  "p(x)  is  that  string 


for  which  the  weighted  average  M Is  a minimum. 


J 
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3.  Structure  of  the  ideal  string. 

Jn  order  to  find  the  density  ^(x)  of  the  ideal  string  we  must  solve  the 
following  rather  atypical  problem  in  the  calculus  of  variations. 

1 


Minimize  M[p] 


1 f(x)  p(x)dx 
^ 0 


(3.1) 


subject  to  the  equality  constraints 

(n  - 1,  2,  . 


n n 


U„[pl  - 
n n 


and  to  the  inequality  constraint 
p(x)  > 0. 


N) 


(3.2) 


(3.3) 


In  order  to  transform  this  optimization  problem  into  a more  standard  form, 
we  shall  first  prove  that  the  ideal  string  must  be  made  up  of  a finite  number  J 
of  concentrated  point  masses,  i.e.  that 


p(x)  •I  m,  6(x-x,) 

j-1  J J 


The  proof  is  of  the  reduct io  ad  absurdum  typ*.  Namely,  if  denotes  the  subset 
of  Ru  which  is  made  up  of  all  the  strings  with  an  infinite  number  of  eigenfre- 
quencies  then  we  shall  assume  that  the  ideal  string  is  in  that  subset,  i.e. 


p(x)eR, 


and  then  see  that  this  assumption  least  to  a contradiction. 
We  start  by  writing 


(3.5) 


p(x)  - [r(x)]‘ 


(3.6) 


a2. 


thus  satisfying  the  inequality  constraint  (3.3).  Since  r (x)eR^  , we  can  solve 

A 00 

(2.4)  and  (2.5)  and  construct  the  set  of  functions  (Wj^(x)}j  defined  as  follows: 


^Mote  that  this  class  of  strings  is  not  equivalent  to  the  class  pel, (0,1). 
Indeed,  even  though  p*l-Hnd(x-ls)^Lf (0, 1),  it  has  an  infinite  number  of  eigenfrequencies, 
A rigorous  treatment  of  this  point  would  require  the  use  of  Stieltjes  integrals. 


± 
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W2n.,(x)  - [0„(x)]‘ 


W2n(x) 


A ^ 

[u„(x)]2 


1,  2, 


(3.7) 


This  set  of  squares  of  the  eigenfunctions  is  complete  (Borg  19^9»  p.  61,  Levitan 
1952, I96A),  This  resuit  is  closely  related  to  the  theorem  regarding  the  uniqueness 
of  the  solution  of  the  inverse  Sturm-Liouvi 1 le  problem,  Heuristical ly,  this  result 
can  be  understood  as  follows.  Consider  two  strings  p(x)  and  p(x)  + 6p(x), 
where  6p  is  a small  density  variation.  The  corresponding  eigenvalues  differ  by 

f 1 


6\  - 

n 


and 


(5p(x)  u (x)dx 
0 n 


p(x)  u^(x)dx 


5p(x)  v^(x)dx 


(3.8a) 


i5y. 


(3.8b) 


1 


p(x)  vMx)dx 
n 


respectively.  If  it  were  possible  to  find  a 6p(x)^0,  which  is  orthogonal  to  all 

the  functions  {W^(x)}j'*’,  then  these  two  strings  would  have  the  same  X and  y spectra. 

But  this  is  not  possible  on  account  of  Borg's  theorem  (Borg  19^6).  Hence  the 

functions  {W^(x)}j**,  and  in  particular  {W^(x)}j**,  form  a base. 

Unfortunately,  the  base  thus  formed  is  not  orthonormal.  This  is  a minor 

nuisance  since  a second  base,  bi-orthogonal  to  the  first  is  usually  required.  We 

denote  this  base  by  {n  (x)},”  and  adopt  the  normalization 

n I 


(3.9) 


i1  A A 

W (x)fi  (x)dx  - 6 , 

m n mn 

where  6 is  the  standard  Kronecker  delta.  The  actual  construction  of  the 
mn 

A 

functions  (x)  can  be  carried  out  by  means  of  a Gram-Schmidt  like  procedure. 


^ 2 

Let  us  consider  next  a string  [r(x)  + 6r(x)]^  which  is  in  i.e.  a 
nearly  ideal  string.  As  can  be  seen  from  (3.8),  the  fact  that  this  string  is 

A A 

in  Ry  implies  that  r6r  is  orthogonal  to  the  functions  W (x)  for  n ■ 1 , 2,  . . . , 2N 
N n 

Consequently, 

00  ^ 

r6r  * y a n (x)  (3.10) 

n-2N+l  " " 


where  the  coefficients  are  arbitrary.  Consequently,  the  variation  in  the 


weighted  mean,  viz. 

(5M  - 2 


j: 


f(x)  r 6r  dx  , 


(3.11) 


can  be  written  thus: 

00 

(5M  - 2 y a 
n-2N+1  '■ 


j; 


f(x)fi  (x)dx 


(3.12) 


Now,  since  p(x)  is  the  ideal  string,  this  variation  must  vanish  for  all 


This  requires  that 


f(x)n  (x)dx  *0  n ■ 2N  + 1 , . . . 

•o  " 


(3.13) 


and  therefore  the  function  f (x)  admits  the  following  finite  series  representation 

2N  A 

f(x)  - I fn  W^(x)  . (3.1^) 

n-1 

By  means  of  simple  manipulations,  we  can  also  write 

2N  ^ A 

f (x)cosY+f’ (x)sinY  ■ I f [W  (x)cosy+W  (x)sinY]  . (3.15) 

**  1 n n n 

n»l 

Now,  recalling  condition  (2.9)  we  can  see  that  near  x ■ 1,  the  left  hand  side  of 

the  above  equation  is  a function  nearly  equal  to  1,  whereas  on  account  of  the 

boundary  condition  at  x ■ 1,  the  right  hand  side  represents  a function  nearly 

equal  to  0.  More  specifically,  we  can  always  find  an  interval,  say  (1-e,l), 

over  which  (3.15)  is  false.  Vte  have  reached  a contradiction  which  implies  that 
A 00  A 

p(x)  \s  not  In  • Consequently  p(x)  must  have  a finite  number  of  degrees  of 
freedom  and  hence  be  of  the  form  {3.k). 
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By  making  use  of  this  knowledge  about  the  structure  of  the  ideal  string, 
we  can  transform  the  original  problem  (3.1)  ■ (3-3)  into  a programming  problem. 
Indeed,  by  replacing  (3.^)  into  {3-1)  “ (3-3),  we  can  state  the  new  problem  as 
f o 1 1 ows : 


J 

Minimize  M * ^ f(x.)  m.  , (3.15) 

j-1  J ■> 


subject  to  the  equality  constraints 

G = det(G..m.  ■ 5 . . ) * 0 

n ij  J U 


r = det  (r,.m,  - —6,.)  * 0 , 
n ij  J IJ 


1,  2,  ...,  N (3.16) 


and  the  inequal  it/  constraints 

0 - Xq  < X,  < . . • < Xj^,  - 1 , 

m j > 0 , j ■ 1 , 2 , . . . , J 


(3.17) 


In  (3.16) 


G. 5 ■ G (xj , X. ) 


IJ 


r. . ■ r(x. , X. ) 

IJ  I J 


and  "det”  stands  for  determinant.  Even  though  it  is  not  possible  to  write  the 
solution  of  (3.15)  “ (3.17)  for  a general  function  f(x)  and  general  spectral  data, 
much  is  known  about  such  problems  (see  e.g.  Gass  1969). 


I 
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4.  A special  case;  the  Krein  string 

If  the  function  f(x)  is  chosen  as  follows: 

f(x)  » ^ (xcosg  + sinB)^  (4.1) 

where  c is  determined  by  means  of  (2.9),  viz. 

c=  (cosB+sing) (cos6cosY+sin6cosY+2cos6sinY)  , (^.2) 

then  the  convex  programming  problem  (3.15)  “ (3»17)  can  be  solved  very  easily. 
Note  that  for  this  case,  the  ideal  string  minimize  a linear  combination  of  the 
mass  Mq  and  the  first  and  second  moments  and  M2.  Indeed, 

M[p]  * ^ (M2Cos^B+2MjSinBcos6+MQsin^6)  . (^.3) 

It  is  possible  to  express  this  combination  of  M^,  Mj  and  M2  in  terms  of 
the  given  eigenvalues,  viz. 


2 2 
M2COS  B+2MjSinBcosB+MQsln  g 


sinlS-c)  cosBcosrsinjBn)  . 
cosacosY+s 1 n (a+y) 


(A.i*) 


where  the  prime  indicates  that  the  term  k * n is  omitted.  This  formula,  which 
is  derived  in  the  appendix,  generalizes  the  formula  for  M^  first  given  by  Krein 
(1951.  1952).^ 

Returning  to  the  ideal  string,  we  can  show  that  the  number  of  degrees 

of  freedom  of  this  string  is  equal  to  N.  Indeed,  if  it  were  made  up  of  J>N 

point  masses,  then  It  would  have  J eigenfrequencies  {X  and  {p  } and 

n I n . 


tKrein's  formula  is  obtained  by  setting  a - y “ 0 and  g ■*  tt/2. 
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J 

M [ ^ m.  6 (x-x . ) ] 

j-1  J J 


1 


where  Q is  a constant  related  to  ct,  B and  y.  Recalling  that  the  y and  A eigen- 
values interlace,  it  is  easy  to  see  that  all  the  terms  in  the  above  sum  are 
positive.  Consequently, 


where 


J N 

M[  I m.  6(x-x.)]  > Q [ 
j-’  J n-1 


1 


k-N+l 


(‘>.6) 


On  account  of  the  ordering  of  the  eigenvaiues  y and  X 

n 

p > 1 . 
n 


it  is  obvious  that 


As  a result 


J 

M[  ]]  m.  6 (x-x.)]  > 

j-1  J J 


M[  I m.  6(x-x.')]  . 
•_1  ' * 


(if. 8) 


N N 

In  other  words,  given  the  truncated  spectra  {X  },  and  {y  the  minimum  of 

n I n I 

2 2 

M2  cos  B + 2Mj  sinBcosB  + M^  sin  B is  reached  for  a string  with  N-degrees  of 
freedom.  Since  there  is  only  one  such  string  in  Rj^,  this  ideal  string  is 
uniquely  determined.  For  the  case  a - y “ 0*  8 - 'n’/2  this  string  is  none  other 


than  the  Krein  string  obtained  by  writing  the  rational  fraction 
(1-  as  a Stieitjes  continued  fraction,  namely 


TT  " • 


n 


N 


where  2.j  ■ . 


i 


K 


12  - 


Acknowl edqments 

I am  grateful  to  Professor  R.  L.  Parker  for  several  discussions  which 
have  stimulated  my  interest  in  the  theory  of  ideal  bodies  and  to  Professor 
J.  Bona  for  constructive  criticisms  of  the  mathematical  proofs.  This  research 
was  supported  by  the  Office  of  Naval  Research  under  Contract  NOOOl l*-76-C-003^ 
and  by  the  National  Aeronautics  and  Space  Administration  under  Grant  NSG-727^. 
Reproduction  in  whole  or  in  part  is  permitted  for  any  purpose  of  the  United 


States  Government. 


- 13  - 


References 

Borg,  G.,  191*6.  Eine  Umkehrung  der  Sturm-Liouvi Ileschen  Ei genwertaufgabe, 

Acta  Math..  2§»  1"96. 

Borg,  G.,  191*9.  On  the  completeness  of  some  sets  of  functions,  Acta  Math..  8l,26l*-283. 
Busse,  F.  H.,  1969.  On  Howard's  upper  bound  for  heat  transport  by  turbulent 
convection,  J.  Fluid  Mech..  37,  l*57"l*77. 

Gass,  S.  I.,  1969.  Linear  Programming;  methods  and  applications.  McGraw-Hill, 

New  York  pp.  358. 

Howard,  L.  N.,  I963.  Heat  transport  by  turbulent  convection,  J.  Fluid  Mech., 

W,  405-1*32. 

Howard,  L.  N.,  1972.  Bounds  on  flow  quantities,  Ann.  Rev.  Fluid  Mech..  4, 

473-494. 

Krein,  M.  G.,  1951.  Determination  of  the  density  of  an  inhomogeneous  symmetric 
string  from  its  frequency  spectrum,  Dokl . Akad.  Nauk  SSSR,  76,  345-348. 

Krein,  M.  G.,  1952.  On  inverse  problems  for  inhomogeneous  strings,  Dokl . Akad. 

Nauk  SSSR.  82,  669-672. 

Levitan,  B.  M.,  1952.  On  the  completeness  of  squares  of  characteristic  functions, 

Dokl.  Akad.  Nauk  SSSR.  83,  349-352. 

Levitan,  B.  M, , 1964.  Generalized  Translation  Operators  and  some  of  their 

appl Ications.  Israel  Program  for  Scientific  Translations,  Jerusalem  pp.200. 
Parker,  R.  L.,  1974.  Best  bounds  on  density  and  depth  from  gravity  data.  Geophysics. 

29,  644-649. 

Parker,  R.  L.,  1975.  The  theory  of  Ideal  bodies  for  gravity  Interpretation, 

Geophys.  J.  R.  Astr.  Soc. . 315-334. 

Sabatier,  P.  C.,  1977a.  Positivity  constraints  in  linear  inverse  problems  - 1. 

General  theory,  Geophys.  J.  R.  Astr.  Soc..  48,  415-442. 

Sabatier,  P.  C.,  1977b.  Positivity  constraints  in  linear  inverse  problems  - II. 
Applications,  Geophys.  J.  R.  Astr.  Soc..  48,  443-460. 


Cambridge. 


r 


t 


- - 


Ti tchmar 
Whi ttake 


1962.  Eigenfunction  expansions;  Part  I.  203  pp.  Oxford, 
Watson,  G.  N.,  1952.  A course  of  modern  analysis,  6O8  pp. 


Appendix 


A generalization  of  Krein's  formula. 

We  shall  be  concerned  with  a string  of  density  p(x)  vibrating  in  two 
different  configurations.  The  eigenvalue  problem  for  the  first  configuration 
i s : 

1 1 

u + A pu  *0 
n n n 

u„(0)cosot  - u ’(O)sina  - 0 (A.l] 

n n 

u ( 1 ) COSY  LI  ’ ( 1 ) s i ny  • 0 


whereas  the  second  is 
1 1 


V 

n 


n 


0 


V ’(O)cosB  ■ V *(0)sin6  ■ 0 

n n 


(A. 2] 


v^ ( 1 ) cosy  + ’ ( 1 ) s i ny  ■ 0 

Just  as  in  the  body  of  the  paper,  we  can  assume  without  loss  of  generality 
that  0<a<6<-y.  As  a result,  system  (A.l)  is  stiffer  than  system  (A. 2) 
and 

U,  < X,  < . . . < < . . . (A. 3) 

In  order  to  solve  (A.l)  and  (A. 2),  we  introduce  two  fundamental  solutions 
yj  (x;X)  and  y2(x;X)  of  the  equation 

II 

y + Xpy  *0  (A. 

such  that 


y,  (1;X)  - -siny  , 
y,’(l;X)  - cosy  , 

y2(’;^)  ■ -cosy  , 
y,’  (1 ;X)  - -siny  . 


(A.  5] 


(A.  6) 


In  view  of  the  linear  independence  of  these  solutions  and  of  the 
conditions  (A. 5),  (A. 6)  we  can  write 

y,(x;X)  y^’(x;X)  » 1 . (A. 7) 

It  is  convenient  to  introduce  two  pairs  of  auxiliary  variables,  namely 
U(x;X)  * cosa  yj(x;X)  - sina  yj’(x;X) 

(A.  8) 

V(x;X)  » cosB  yj (x;X)  - sinB  y,’(x;X) 

and 

$(x;X)  - cosa  y2(x;X)  - sina  y2’(x;X) 

(A. 9) 

ffxjX)  - cos3  y2(x;X)  - sing  y2’(x;X) 


Several  remarks  are  now  in  order.  First  of  all,  for  x * 0 the  zeros  of 
U and  V are  the  eigenvalues  of  (A.l)  and  (A. 2).  Consequently,  in  view  of  the 
fact  that  U(0,X)  6 V(0,X)  are  entire  functions  of  X of  order  1/2  (Titchmarsh 
1 962) , we  can  wr i te 


U(0,X) 


u(o,o)  yj  (1 

n“I 


(A. 10) 


V(0,X)  - V(0,0)  TT  0 " “ ) • 

n«l 

Next,  we  should  point  out  that  the  Wronskian  equality  (A. 7)  can  be  written  in 
terms  of  the  new  variables  as  follows: 

U(x;X)  - V(x;X)  $(x;X)  ■ sin(a-6)  . (A. 11) 

Let  us  now  consider  the  function 


F(z-i)  . v(0;0)  non)  - Ho-,0)  v(o;2) , 

’ ' z(z-x)  v(o:Z) 

This  is  a merotnorphic  function  with  simple  poles  at  Z ■ X and  Z ■ (n-1 ,2, . . . ) but 

not  at  Z - 0.  If  is  a circle  in  the  Z-plane  of  radius  |Z|  - X^,  then  It  is 


possible  to  show  that 


-A. 3 - 


Several  simplifications  are  possible.  In  particular,  since  V(0;y^)  * 0,  we  can 

exploit  (A. 11)  to  write  '?(0;y  ) in  terms  of  U(0;y^);  also  we  can  use  (A. 10)  to 

n n 

evaluate  3V(0,X)/3X.  Therefore 


So  far,  we  have  just  repeated  the  various  steps  in  the  proof  of  the  Mi ttag-Leff ler 
theorem  (Whittaker  6 Watson  1952).  We  now  let  X 0.  The  above  formula  becomes: 


U(0;0)  Aim  1 { V(0;0)'1'(0;X)  - ’i'(0;0)V (0;X)  } 

V{0;0)  X-K)  I 


- sin(a-6)  I 
n“l 


(A.  14) 


We  shall  now  see  that  the  left  hand  side  can  be  expressed  in  terms  of  the  mass 
and  first  two  moments  of  the  density  distribution.  To  that  effect,  let  us 
expand  y.(x;X)  in  powers  of  X and  compute  the  first  two  terms  by  substituting 
in  (A.l)  6 (A. 2).  Omitting  the  intermediary  calculations,  we  get 
y^  (0;X)  ■ -cosY“SinY+X[(Mj-M2)cosY+M^sinY]  + . . . 


y^’(0;X)  ■ COSY+  X[-{MQ-M^)cosY”MQsinY]  + . . . 


(A.  15) 


- A.k  - 


and 

■ sinY-cosY+X[- (Mj-M2)sinY+M^cosY]  + . . . 
Y2’(0;X)  * sinY+X[MQ-M^)sinY-MQCOSY]  + . . . 


Replacing  these  expressions  in  (A. 8),  (A. 9)  we  deduce  that 
V(0;X)  ■ -[cosBcosy  + sin(8+Y)] 

+ X[(M^-M2)cos6cosy  + M^cosBsinY 

+ (MQ-M^sinBcosY  + M^sinSsinY]  + . . . 


and 


’i'(0;X)  * [cosBsinY  " cos(6+y)I 

+ A[-(Mj-M2)cos6sinY  + M^cosBcosy 

- (MQ-M^sinSsinY  + M^sinBcosYl  + . . . 


and  consequently 

lim  1 {V(0;0)'1'(0;A)  - ^(0;0) V (0;X) } - (M,-M,)cos^B 

X^O  X ' ^ 

- (cosB+s inB)cosB  + (MQ-M^sinBcosB 

- MQ(cosB+sinB)sin6 

As  a result  (A. 14)  reads: 


cosacosY'*‘S  i n 
cosBcosy+s i n (B+y) 


(M2C0S^B  + 2M^sinBcosB  + M^sin^B)  ■ 


sin(B-a)  I 

n-1 


For  a“0,  6*^/2  and  y^O,  the  formula  reduces  to 


which  is  the  formula  given  by  Krein  (1952). 


(A. 16) 


(A.  17) 


(A. 18) 


(A. 19) 


(A. 20) 


